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ABSTRACT 

The separation distribution for M-dwarf binaries in the AstraLux survey is narrower 
and peaking at smaller separations than the distribution for solar-type binaries. This 
is often interpreted to mean that M-dwarfs constitute a continuous transition from 
brown dwarfs (BDs) to stars. Here a prediction for the M-dwarf separation distribution 
is presented, using a dynamical population synthesis (DPS) model in which “star- 
like” binaries with late-type primaries (< I.5M0) follow universal initial distribution 
functions and are dynamically processed in their birth embedded clusters. A separate 
“BD-like” population has both its own distribution functions for binaries and initial 
mass function (IMF), which overlaps in mass with the IMF for stars. Combining these 
two formation modes results in a peak on top of a wider separation distribution for late 
M-dwarfs consistent with the late AstraLux sample. The DPS separation distribution 
for early M-dwarfs shows no such peak and is in agreement with the M-dwarfs in 
Multiples (MinMS) data. We note that the latter survey is potentially in tension with 
the early AstraLux data. Concluding, the AstraLux and MinMS data are unable 
to unambiguously distinguish whether or not BDs are a continuous extension of the 
stellar IMF. Future observational efforts are needed to fully answer this interesting 
question. The DPS model predicts that binaries outside the sensitivity range of the 
AstraLux survey remain to be detected. For application to future data, we present a 
means to observationally measure the overlap of the putative BD-like branch and the 
stellar branch. We discuss the meaning of universal star formation and distribution 
functions. 

Key words: binaries: general - stars: low-mass - stars: late-type - stars: formation 
- stars: kinematics and dynamics. 


1 INTRODUCTION: COMMON OR 
SEPARATE POPULATIONS? 

Binaries are a do mi nant ch annel of star fo r mation 
jGoodwin_^Krou 2 a l2005l : iDuchene &: Krau^ l2013l : 


iReipurth et al 


20141 1 which is a result of the angular- 
momentum problem in star formation. The distribution 
of their orbital separations are commonly described by 
bell-shaped distributions that depend on the spectral-type 

of the primary. _ _ 

The pioneering work of iDuauennov fc Mayor! ill 99 il l 
showed a wide distribution peaking at about « 30 AU, and 
a binary fraction of 58 ± 10% for G-dwarf binaries in the 
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Galactic field. More recent analyses largely con firmed their 
findings llRaghavan et al.ll2010l : lTokovininll2014l l but suggest 
a some what lower binary fractio n of « 46 ± 2%. Subsequent 
work bv lFischer fc Marcvl (Il992ll suggested a similarly wide 
distribution for the M-dwarf binaries which peaks at a sepa¬ 
ration comparable to solar-type stars. Their binary fra ction 
was estimated to be « 42 ± 9%. iDelfosse et al.l ll2C)04l l’s M- 
dwarf sample confirmed a wide dis tribution. Similar results 
are obtained for K-dwarf binaries llMavor et al.lfl993l . 

Studies of BD binaries later showed the distribution of 
their orbital separations to be significantly narrower and the 
binar y fraction of around ~ 15% to be lo wer than that for 
stars llClose et al ] |2003l: lBouv~ et al. |[20 o 3). This has led to 
the suggestion that BDs are a population separate from the 
hydrogen-burning stars. 
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Due to the different behav iour of stars and BDs in terms 
of their pairing characteristics. pThies fc Kroupal (l2007l . [2()08l f 
suggested an IMF which has a discontinuity around the 
hydrogen-burning mass limit, but with a non-negligible over¬ 
lap of both populations, i.e. some BDs may form “star-like” 
while some very low-mass (VLM) stars form “BD-like”, i.e 
in a fragmenting circumstellar disc. In this model both pop¬ 
ulations pair their objects among each other to form bina¬ 
ries but do not mix in pairing. The two-component IMF 
obtained after correcting the star counts for unseen com¬ 
panio ns is consistent with observations llXhies fc Kroup^ 
I2OO8II and smoothed particle hydrodynamical (SPH) compu¬ 
tations which produce BDs fr om encounter-trigg ered pertu- 
bations in circumstellar discs llThies et al.ll201^ . The mod¬ 
els of Thies fc Kroupa furthermore support an observed 
BD desert (iMcCarthv et al.l I2OO3I : iGrether fc Line weaver! 
l2006l : [Dieterich et al.ll2012l) . Additionally. Thies et al. 1 2015l l 
recently demonstrated that assuming BDs to be a con¬ 
tinuous extension to the hydrogen-burning stars leads to 
an overestimation of the numbers of BDs when com- 
p ared with actual ob s ervat io ns, conh r ming the deductions 
of lKroupa fc Bouvierj ll2003h . [Li et al.l (l20ld ~l recently inves¬ 
tigated the outcome of a similar disc fragmentation mode 
(through self-instability of the disc, though, not through 
encounter-induced fragmentation). By combining SPH and 
Nbody techniques, they found a narrow separation distri¬ 
bution for BD-like objects peaking between 5 and 10 AU 
after 10 Myr of evolution (their fig. 11), similar to the BD 
Galactic held properties, and within the limits of their ini¬ 
tial conditions, their results are in good agreement with the 
observed BD desert. 

The distinguishing characteristics of the two popu¬ 
lations are now being challenged by recent observations 
of M-dwarfs in the Gala ctic held. The AstraLux survey 
lljanson et al ]|20lll2014l') suggests that the M-dwarf separa¬ 
tion distribution is signihcantly narrower than that of solar- 
type sta rs. It is unlikely that a sola r -type separation distri - 
bution llDuauennov fc Mayor! Il99ll : iRaghavan et al.l I 2 OI 0 I I 
is a parent function for both early and late M-dwarfs. 
This is interpreted as evidence for continuous star forma¬ 
tion from BDs to stars, where a population’s binary frac¬ 
tion increases smoothly with increasing mass of the primary 
component, alon g with the peak and wid t h of the separa¬ 
tion d istribution i Duchene fc Kra^l2013fl . IParker fc Meverl 
(l2014ll argue that these observati ons in combination with the 
multiplicity pro perties of BDs ( Thi es fc Kroupal l2007f l. of 
solar-type stars havan et aklboioll and ofXtvpe stars 
dPe Rosa et al ]|2014|) disfavour universal initial binary prop¬ 


erties throughout these spectral types. Instead they consider 
a scenario in which the field binary population is indicative 
of the primordial multiplicity conditions in the star forma¬ 
tion regions they were born in. This would violate the univer¬ 
sality hypothesij^ of star formation according to which the 
universal ity of the IMF and binary populations are st rongly 
coupled dKrou ^l200ll : iKroupa fc Petr-Gotzensl[201lh . 

The influence of a separate BD-like population super¬ 
posed onto a star-like population on the multiplicity prop- 


^ Note that the universality hypothesis of star formation refers 
to an environment independence and does not preclude a mass- 
dependence (see Sec. [ 5 :^ . 


ns 

o 

6) 

o 

< 


Z) 

< 


OS 

o 


0.2 - 


0.15 


0.1 - 


0.05 - 


brown dwarfs (logioapeak=0'66, o=0.4, f[,=15%) 
early AstraLux (logi()apej,|^=1.2, cr=0.8,1^=25%) 
late AstraLux (log,Qap..;,=0.78, cj=0.47, fh=18%) 
MinMS (logi(,apea|^=0.^ a=1.34, L=25%) 
solar-type (logioapgg|j=1.64, a=1.52, f(,=46%) 





H' 



•«' ■ 




! \ 



'4 

’! ■ 

■ 


'Ml 


■ 


1 '1 


■ 

■ 


1 1 

i| 




0 
0.01 




0.1 


1 10 100 1000 10000 

semi-major axis a / AU 

Figure 1. Observed orbital separation distributions and Gaus¬ 
sian fits for binaries in the Galactic field. M-dwarf data from 
the AstraLux survey (histograms and corresponding same- 
linetype Gaussian distributions, statistically cleaned samples 
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2012|.|2014^ and the MinMS data Oower limit. 

IWard-Duone et al. 

201511. BD data from the VLMBA archive as 


IRaghavan et al.l ll201Clh . The key provides the distribution peak 
(logjQ(apeak/AU)), the distribution width (cr) and the binary 
fraction (/u) within the respective sensitivity ranges as stated 
by the authors. 


erties of M-dwarfs has hitherto not been tested. This is the 
aim of the present research paper. In Section [2] the charac¬ 
teristics of the here used surveys are summarized. Section [3] 
outlines how the Galactic field population is modelled. The 
results of our analysis are given in Section |T] and discussed 
in Section [S] The paper closes with a summary in Section!^ 


2 SURVEYS 
2.1 AstraLux 


In a lucky-imaging survey the AstraLux survey in- 
vesti gates binarity amo ng “early” (spectral types MO- 
M6, Janson efcaf l2012h and “late” M-dwarfs (M2-M8, 
Ijanson et al. 2014 1 in the Galactic held. Although the two 
samples are not distinct in mass we nevertheless henceforth 
refer to “early” and “late” for t he sample with the smaller 
l|janson et al.ll2()12 1 and larger djanson et al.ll2()14[l median 
spectral type, respectively. 

Constructing two samples distinct in spectral type or 
primary mass would be desirable in order to amplify the ob¬ 
served mass-dependency under study. The selection criteria 
for the two samples were however systematically different 
(aside from the different spectral type ranges). The early 
sample was selected on the basis of X-ray brightness, while 
the late sample was based on more nearby stars with a cut¬ 
off in infrared brightness. This means that the completeness 
of the surveys are systematically different, and it also follows 
that the early sample is systematically young and often pre- 
main sequence, while the late sample is not. This in turn 
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means that the stellar masses and the spectral types the 
different samples correspond to are systematically different 
as the temperature evolves with age until the main sequence. 
Hence, a cut in spectral type between two sub-populations 
in a combined sample would be inaccurate. The youth of the 
early sample additionally implies that masses of some tar¬ 
gets are still subject to change which might further diminish, 
or amplify, the overall offset in primary mass. Furthermore, 
the majority of stars in the late sample have parallactic dis¬ 
tances while the majority of stars in the early sample do 
not, so adding stars from the early sample to the late sam¬ 
ple would compromise the quality of the semi-major axis 
distribution determination of the latter. This is why we will 
focus on the late-type sample whose separation distribution 
is more reliable. Our conclusions are not affected by inves¬ 
tigating the two samples individually since the DPS model 
(Sec. O takes into account the individual and overlapping 
mas s-ranges. _ 

Ijanson et al.l (120141 ') define a “statistically cleaned” 
(SC) sample. If the combined magnitude of the binary com¬ 
ponents exceeded the limiting survey brightness, the authors 
removed the target from their statistical analysis. Had the 
primary been a single object, it would not have been se¬ 
lected. This SC sample then has 48 binaries among 268 tar¬ 
gets, corresponding to a binary fraction of 17.9% within the 
sensitivity range, logj^g(a/At/) = —0.4 to 1.95, where a is 
the semi-major axis. Fig.[T]shows that, in c omparison to the 
separ ation distribution for solar-type stars jRaghavan et al.l 
I 2 OI 0 I I. both early and late M-dwarf samples exhibit sepa¬ 
ration distributions which are rather narrow around their 
peak. For the two AstraLux samples in Fig. [T]we note the 
following two points: 

(i) The semi-major axis distribution of late M-dwarfs is 
strikingly close to that for BDs. 

(ii) The amplitude of the peak in the late M-dwarf semi¬ 
major axis distribution is larger than the one for early M- 
dwarfs. 

At first glance, item (i) appears to favour the continuity of 
star formation over and above the hydrogen-burning mass 
limit. However, it will be demonstrated that the origin of 
the observed proximity of both distributions might be that 
a separate BD-like population overlaps with a population of 
star-like bodies. In addition, the here devised DPS model 
yields item (ii). 


2.2 MinMS 

The 15 pc volume-limited M-dw arfs in Multiples survey 
(MinMS. IWard-Duong et al.ll2015l ') detects K7-M6 binaries 
over a separation range logjQ(a/j4[/) ~ 0.5 — 4, using an 
infrared adaptive optics (AO) technique to find close com¬ 
panions (~ 10° — 10^ AU) and digitised wide-held archival 
plates for wide companions (10^ — 10^ AU) each covering 
multiple epoch s. The targets stem from a reduced Hippar- 
cos catalogue Jvan LeeuwerJl200^ and have reliable paral¬ 
laxes. 65 co-moving stellar companions in a sample of 245 
late-K to mid-M dwarfs yield a companion star fraction of 
23.5 ± 3.2 per cent within the sensitivity range. 

The frequency of detected binaries with a given separa¬ 
tion increases towards smaller separations (Fig. [T]). As mea¬ 
sured by the Gaussian width, the ht to the MinMS separa- 



Figure 2. The stellar and substellar IMF are disjoint (black his¬ 
tograms). The dotted lines indicate an underlying power-law of 
the form with a = 0.3, 1.3 and 2.3 (from left to right). “BD- 
like” VLM stars are objects from the BD branch of the IMF lying 
above the hydrogen burning mass-limit (vertical dashed line) and 
“star-like” BDs stem from the stellar branch of the IMF lying be¬ 
low the hydrogen burning mass-limit. Both populations mix and 
form the system IMF (31% and ^ 100% binary fraction for the 
BD and stellar branch, respectively, see Tab.[T]l. The discontinuity 
is thus hidden from an observer (grey histogram). 


tion distribution in IWard-Duong et al.l ll2015h suggests that 
it is at least 1.7x wider than for the early AstraLux separa¬ 
tion distribution, although the covered spectral-type ranges 
are comparable in both surveys. This potential tension will 
be qualitatively addressed in the discussion of our results 
(Sec. HI . 


3 A MODEL FOR THE GALACTIC FIELD 


The AstraL ux M-dwarf dja nson e t al.l I2OI2I . l2014l ) and 
MinMS data llWard-Duong et al.ll2015l ) are for Galactic field 
binaries. A model is thus needed to describe how these em¬ 
anate from their birth embedd ed stellar clusters. Such a DPS 
model is readily available from iMarks fc Krounal (|201ll . Sec- 
tion l3.3l) . A separate BD-like population fSection l3.2l) is here 
added to a star-like population. For the star-like population, 
we assume universal initial conditions for late-type stars, 
and then account for the subsequent dynamical processing 
(Section [3T} . 

The DPS model will provide predictions for M-dwarf 
binaries in the Galactic field. These predictions do not stem 
from first principles because a theory of star formation does 
not exist which yields such information. Instead, the bi¬ 
nary orbital-parameter distributions to start with before dy - 
namical processing result from inverse DPS (lKroupalll995h 
of the then available emp irical data for solar-type stars 
(iDuquen nov fc Mavoilll99ll) and pre-main sequence binaries 
(e.g. Leinert et al.lll993 ). But since the involved DPS model 
parameters (Tab.[T]below) have been constrained using inde¬ 
pendent samples, as explained in the forthcoming sections, 
the DPS model with a set of previously constrained pa- 
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rameters can be tested against other observations using the 
same parameters. This is done here for the AstraLux and 
MinMS data. 


3.1 Universal formation plus dynamical 
processing for late “star-like” objects 

As a working hypothesis to address the question of the 
(non-)continuity of star formation the birth binary popu¬ 
lation for star-like objects in embedded star clusters is as¬ 
sumed to be the same for all late-type stars (iKrou ^ l 201 ll: 
iKroupa &: Petr-Gotzenj|201ll l. “Star-like” refers to all bod¬ 
ies selected from the two-part IMF which extends into the 
BD regime (producing star-like BDs down to O.OGM©) and 
overlaps with a separate “BD-like” populat ion (see Fig. [2] 
and Section iThies fc Kroupall20CI7l . l2008l ~). 

The birth binary population is built from an universal 
separa tion distribution for late-type stars derived bv lKrou^ 
ll 199511 with an excess of long p eriod (soft) b inaries over 
that in the Galactic field (eq. 8 in iKrou ^ Il995h , as seen in 
young star formation regions such as Ta urus Auriga (Tau; 
iKohler fc Leineri 19981: iKraus et al.ll201lll and Corona Aus- 
tralis (CrA. iKohler et al.ll2008l l. The initial binary popula¬ 
tion is defined using random pairing of binary component 
masses for low-mass stars (<a few Mq), which are selected 
from a standard two-part power-law stellar IMF llKroupal 
I2OOI), along with a th ermal eccentricity distribution (see 
also Kroupa et al.ll2013ll which bends ov er in the course of 
dynamical processing ( Marks et al their figs. 2 and 5). 

While random pairing does not well reproduce the rather flat 
mass-ratio distribution for solar-type stars and more massive 
primaries llReggiani fc Mev^l201lf : iDuchene fc Krau^l2013l l 
it is in agreement with the weakly rising q-distribution for 
M-dwarf binaries and also with the over all mass ratio dis¬ 
tribu tion of late-type stars (see fig. 5 in iMarks fc Kroup^ 
[20 i 3). The results on the semi-major axis distribution pre¬ 
sented here are not affected by the chosen pairing mecha¬ 
nism since dynamical binary processing shows only a weak 
preference to disrupt low -q systems first llMarks et al.ll20jT] : 
IParker fc Reggianill2013l '). Any additional orbital parameters 
are calculated using Kepler’s laws. The subsequent birth bi¬ 
nary p opulation is su bjected to pre-mainsequence eigenevo- 
lution llKroupal Il995h to account for gas-induced changes 
of orbital parameters in the c ircum-binary material during 
the c luster formation (see also lStahlerj|201(il : iKorntreff et all 
I2OI2I I. The initial population, which is prone to dynamical 
processing inside its host environment, i.e. an embedded star 
cluster, is thus arrived at. 

In order to calculate the effects of dynamical process¬ 
ing on initially binary-dominated populations inside young 
star clusters, we apply the analytical description of the 
dynamical p r ocessi ng in Nbody computations devised by 
iMarks et al.l (l201ll l. In terms of the resulting binary pop¬ 
ulation, the outcome should be dynamically equivalent to 
an initi ally sub-struc t ured and dynamically cold configura¬ 
tion, as IParker et al.l (l201lll conclude. The initially binary- 
dominated population changes on a crossing-time scale, i.e 
the resulting distribution aft er a few Myr depe nds only on 
the initial stellar density fsee IMarks et al.ll 201 lll . 


3.2 A separate “BD-like” population 

Given the existence of the BD desert and the apparently 
different binary characteri st ics be t ween stars and BDs (Sec¬ 
tion [1]), iThies fc Kroupal (I2OO7I . I2OO8II quantified a stel¬ 
lar plus BD IMF which has a discontinuity around the 
hydrogen-burning mass limit. However, the stellar and BD 
parts have an overlap (Fig.[2|, since the star formation pro¬ 
cess does not care about hydrogen burning and it is unclear 
why the hydrogen-burning mass limit should constitute a 
sharp cutoff. Note that upon combining both populations an 
observer would not readily see the discontinuity but rather 
a declining and cont inuous IMF around th e star-BD mass 
limit (see Fig. l^and lThies fc Kroupall2007ll . 

In this contribution a “BD-like” population is added 
to the star-lik e population result i ng fro m the Galactic field 
DPS model of IMarks fc Kroupal (1201 ll . Section [3.311 . A BD 
or late M-dwarf is called “BD-like” if it stems from the BD- 
like branch of the IMF (masses between 0.01 and 0.2Mq, 
Fig. 121). The exact value of the upper cut-off mass for BD- 
likes is not k nown and might vary between r egions. The best- 
fit values in iThies fc Kroupal 1I2OO7I . I2OO8I I for Tau, IC348 
and the Trapezium cluster in the ONC are in the range 
0.1 — 0.2Mq, with individual uncertainties of about ±0.04 
to ±0.08 (their table 3). For Tau its value might be closer 
to O.IM 0 , and to 0.2Mq for IC348 and in the ONC, but 
their values agree within their 2a mutual errorbars. Since 
Tau type aggregates cannot be dominant in contributing to 
the Galactic field, give n a significant soft binary component 
not p resent in the field dMarks fc Kroup^l 2012 l : TMarks et al.l 
I20141 . we here adhere to rrimax.BD = 0.2Mq determined for 
IC348 whose binary populatio n resembles the Galacti c field 
at long periods more closely dMarks fc Kroup3l2012fl . The 
most massive objects forming in th e discs around host stars 
of 0.7Mq used in the analysis of iLi et al.l d2015lf is trun¬ 
cated at « 0.2Mq as well, motivated by the mass spectrum 
of objec ts emerging in circumstella r discs in SPH computa¬ 
tions of IStamatellos fc WhitworthI d2009fl . The influence of 
the choice of the BD-like cut-off mass is discussed further in 
Sec. 15.11 Note that our chosen value must not be taken rep¬ 
resentative for any star forming region since we discuss here 
a BD population originating and superposed from many star 
forming events. 

The BD-like population is excluded from dynamical 
processing in its birth cluster (as is described in Section [3.II) . 
since we do not have good constraints on the initial BD 
population at this time. Note that this does not imply that 
we consider this population to be dynamically inactive, i.e. 
we do not assume the field population to resemble the pri¬ 
mordial population. Instead the dynamical processing of the 
BD-like part is implicitely accounted for by constructing it 
such that upon superposing it with the dynamically pro¬ 
cessed star-like population, the BD binaries (< O.O 8 M 0 ) 
added from both branches match the BD binary charac¬ 
teristics in the field. This is possible as the outcome of any 
dynamical model has to reproduce what is observed for BDs 
in the Galactic field as a constraint. This approach, however, 
fixes the parameters of the BD-like population and the dy¬ 
namical population synthesis parameters a priori (Tab. [T]). 
The here used BD-like population is thus the one for the 
Galactic field and must not be used for population synthesis 
in star clusters. The BD-like population that should act as 
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the input to star cluster models could be universal, but here 
we do not address this issue. 

To model BD-like binaries, individual masses are se¬ 
lected from the BD IMF (Fig. [5]) randomly and stored in 
an array. From this array, masses are paired such that the 
observed mass-ratio distribution for BDs in the field, which 
has a strong preference toward unit mass-ratio, is repro- 
duced. To d o so, t he biased-pairing algorithm introduced by 
iThies et al.l (l2015l i is used. This procedure applies a prob¬ 
ability p — to each binary, where q = m 2 /mi £[ 0 : 1 ] 
is the mass ratio and 7 = 2. If another random number 
is smaller than p the binary is accepted. If it is larger a 
new companion for the first selected object is assigned and 
the procedure is repeated. This ensures that a mass-ratio 
close to unity is strongly preferred for BDs. Note that a re¬ 
jected companion is not discarded. It later becomes either a 
component of a different binary or stays single. This step is 
important to maintain the shape of the BD-like branch in 
the single-star IMF (Fig. [2)). Note also that biased pairing is 
a natural outcome of t he formation of B Ds in fragmenting 
stellar accretion discs llThies et al.ll201^ . Each BD-like bi¬ 
nary is given a semi-major axis selected from the observed 
separation distribution of BDs in the Galactic field 0 Most 
of the selected BD-like objects remain single to yield the 
observationally constrained binary fraction of « 15 per cent 
for BDs upon combining it with the processed star-like pop¬ 
ulation. The total binary fraction for BD-like binaries in the 
Galactic field is then ~ 31%. The discontinuity in the IMF 
is measured by the fraction of BD-like to star-like objects 
in the population, TZpop = N'BD-iike/N'star-iike and chosen 
to be 0.3. This reflects the empi rically determined valu e for 
Taurus-Auriga and the Pleiades (iThies fc Kroup^l2007l l. i.e. 
for every three star-like bodies there is one BD-like object. 

3.3 Galactic field dynamical population synthesis 

The stellar and BD population in a galaxy is the result 
of the addition of all populations formed in all embed- 
de d clusters. According to the galactic field DPS model 
of [Marks fc Kroupal (l201ll '). each embedded star cluster’s 
binary population is processed for 3 Myr, the time at 
which the first supernovae are expected to occur and in 
so doing drive out the resid ual-gas and destroy most of 
their natal stellar aggregate (iLada fc Ladal [ 2 OO 3 I 1 FI These 
clusters are thus the building blocks of the stellar sin¬ 
gle and binary population of the Galactic field. Embed¬ 
ded stellar cluster masses, M^d, are assumed to be dis¬ 
tributed according to a single power-law initial embed¬ 
ded clus ter mass function (ECMF), oc Mfif , with 
/3 ~ 2 jde la Fuente Marcos fc de la Fuente Marcod |2004 

iGieles et al.ll200^ . The dynamically processed stellar popu¬ 

lation of a cluster with a given initial (embedded stellar) 
mass is weighted (multiplied) by the number of clusters 
that have this mass according to the ECMF. The sum of 
all these weighted popul ations results in the Ga lactic field 
population. As shown by [Marks fc Kroupal ll201lll . the DPS 

^ See Section ED for further discussion of this assumption. 

® The results are not very sensitive to this time-span as dynamical 
processing of initially binary dominated objec ts occurs rapidly o n 
a crossing-time scale of the embedded cluster llMarks et al.l20llh . 


Table 1. Parameters defining the initial star-like population in 
embedded clusters, the BD-like population in the Galactic field 
and the Dynamical Population Synthesis procedure. From top to 
bottom the parameters denote (for details see the text): initial 
binary fraction and minimum mass for objects on the star-like 
branch, field binary fraction and maximum mass for objects from 
the BD-like branch, fraction of BD-like to star-like objects and 
biased pairing exponent, star formation rate in the Milky Way, 
power-law index of the ECMF, minimum embedded star cluster 
mass and typical embedded cluster half-mass radius. 


Parameter 

Value 

star-like population (embedded clusters) 

fb 

Ri 100% 

^min 

0.06 Mq 

BD-like population (Galactic field) 

fb 

31% 

^max 

0.2 Mq 

^pop 

0.3 

7 

2 

Dynamical Population Synthesis 

SFR 

3 Mq yr ^ 


2.0 


5 Mq 

rh 

0.1 pc 


model works well to simultaneously describe various orbital- 
parameter distributions of the Galactic field binary popula¬ 
tion with a single set of parameters, depending on spectral- 
type (mass) of the primary component. This notion is sup¬ 
ported by the present work and will be further substantiated 
in an upcoming contribution, which uses the latest observa- 
tional data. In this wor k the DPS parameters inferred by 
[Marks fc Kroupal |[201ljl will be used (Tab.[T|). 


4 RESULTS 

4.1 Separation distribution 

Using a star-like formation mode only does not match the 
recent AstraLux data for late M-dwarfs. Fig.[3](left panel) 
however demonstrates that upon adding BD-like M-dwarfs, 
i.e. those late M-dwarfs which form as part of the BD branch 
in Fig. (2] a peak on-top of the wider star-like distribution ap¬ 
pears. This is due to, in the present formulation, the BD-like 
M-dwarfs sharing their separation distribution with those for 
BDs in the Galactic field. The DPS model binary fraction 
of 21.5% in the range logj^Q a = —0.4 to 1.95 is in reasonable 
agreement with the 17.9% observed in the late M-dwarf SG 
sample. 

For a more sophisticated analysis, we compare the dis¬ 
tribution shapes by means of a Kolmogorov-Smirnov (KS- 
)test. We do this by subjecting the DPS model to the com¬ 
pleteness of the survey. One thousand random realizations of 
the DPS model are generated and the median of the match 
probability is adopt ed. This is the same method used by 
[.lanson et al.[ l[2014l l to find a Gaussian distribution that 
best describes the data. For details see their sec. 6.2. The 
test involves the DPS model mass-ratio distribution which is 
discussed in Section [121 The experiment results in a match 
probability of 30.4%, i.e. the hypothesis that the DPS model 
constitutes a parent distribution of the data is confirmed to 
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Figure 3. Comparing orbital semi-major axis distributions of the DPS model (curves) with AstraLux and MinMS data, respectively 
(histograms). For reference, the semi-major axis distribution of solar-type stars (dotted blue line) is shown in both panels. The DPS 
model binary mass ranges are chosen according to the mass estimates of the binaries in the AstraLux and MinMS survey. Line-types 
as in Fig. [T] Left panel: Upon combining the star-like M-dwarfs (dashed black curve) with the BD-like M-dwarfs (0.08 — 0.2 Mq on the 
BD-like branch) in the DPS model, whose distribution resembles the one for (star-like -|- BD-like) model BDs (solid black curve), the 
late M-dwarf DPS model results (green dashed curve). It resembles the late -type data in the observ ed range (green dashed histogram, 
SC sample). For comparison, the M-dwarf data from the pioneering study bv iFischer fc Marcvl lll992l . crosses) is shown. Right panel: A 
peak is not expected (solid yellow and dashed-dotted red curve) in the early M-dwarf data (corresponding histograms) since the BD-like 
branch here extends to O. 2 M 0 only. 


within « Icr certainty. This result is discussed further in 
Section [ 5 l 

It is i nteresting to note that th e now more than 20-year 
old data of iFischer fc Marcvl ll 19921 ) shows a similar excess of 
M-dwarf binaries which is located at about the peak of the 
field BD and late M-dwarf separation distribution. Primary 
masses in their observations were in the range ~ 0.1— 0.6Mq, 
i.e. a contribution of BD-like M-dwarfs is expected. The huge 
uncertainties did not suggest a real feature. But, retrospec¬ 
tively, as a tracer of a separate BD-like population, the peak 
might have been in front of our eyes all along. 

For early M-dwarfs a similar peak on-top of a wider dis¬ 
tribution is not expected (Fig. [31 right panel). This is due to 
the BD-like branch extending, in the present formulation, up 
to 0.2Mq only (F ig. (21). This va l ue co mpares to the lowest- 
mass M-dwarfs in jjanson et al.l (l2014l )’s early sample. Note 
how the different height of the observed AstraLux early 
M-dwarf distribution is reasonably matched by the DPS 
model. The wings of the fitted distributions as seen in Fig. [T] 
are not reproduced, but with the exception of the lowest 
separation bin the DPS model might even compare with 
the raw AstraLux data. On the other hand, the available 
MinMS data and the distribution fitted to this data in Fig.[T] 
is in excellent agreement with the DPS modelQ 

This immediately questions how the observational data 
on early M-dwarfs, which share similar properties (Sec. [2]), 


^ Note that, although MinMS targets with < 0.2 Mq exist, we 
have run the DPS model without a BD population since only one 
binary with < 0.2 Mq is part of the observed separation distribu¬ 
tion. 


relate to one another. While the increased frequency of bi¬ 
naries with distant companions can be readily explained 
through the enhanced coverage of semi-major axes in the 
MinMS sample, the conflict remains for the closest cov¬ 
ered separations (Fig. [T]). This is at the very least puz¬ 
zling since the covered spectral-type range is comparable 
in both surveys. If the slightly wider mass-range of the 
MinMS targets is responsible, yet unknown processes could 
be at work around the M-/K-dwarf boundary. Another pos¬ 
sible caveat is the systematic age difference between the sam¬ 
ples (Sec. l2.lT) . However, at present it is hard to see a reason 
how and why these differences should approximately double 
the observed distribution width. 

Alternatively, one may ask whether the fits of Gaussian 
distributions to the AstraLux and MinMS data are reli¬ 
able representations of the underlying parent distributions. 
The raw separation distributions (histograms in Fig.O right 
panel) appear to compare better in the separation range 
where both studys overlap than a comparison of the fit¬ 
ted distributions suggests. Is it then possible that both dis¬ 
tributions are not distinct after all and stem from a com¬ 
mon parent distribution? In this latter case the Gaussian 
fits would be too simplistic. We cannot decide this on the 
basis of the present study and pos sibly requires f urthe r ob- 
servations. It is noted tha t neither jjanson et al.l (|2012| ) nor 
IWard-Duong et Id] (l2015l ) attempt to make any detailed fit 
involving thorough statistical tests but merely discuss var¬ 
ious general options. So it is as of yet not clear whether a 
parent distribution exists that fits both datasets simultane¬ 
ously at acceptable confidence limits. A wider distribution 
for early M-dwarfs as in the MinMS da ta compares bet- 
ter with some previously obtained results (iFischer fc Marc^ 
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Figure 4. Oberserved and DPS model mass-ratio distributions 
for binaries in the semi-major axis range log(a/AC/) = —0.4 
to 1.95, as in the observations. Combining star-like (dashed grey 
line) and BD-like M-dwarfs (solid grey line) as in Fig. [3] the 
late M-dwarf DPS model appears (solid black curve). The his¬ 
togram with Poissonian errorbars shows the late AstraLux data 
(SC sample). Caution is required since mass estimates for 
AstraLux targets are uncertain fSec. [4?^ . 


I 1992 I : iDelfosse et al.l |2004| ) , while a narrower distribution 
as in the early Ast raLux data is suggested as well by 
iBergfors et alJ (l2niflh ’s study. Note that the latter two stud¬ 
ies which find a narrow distribution used the same observa¬ 
tional technique, i.e. lucky-imaging. 


4.2 Mass ratio distribution 

Masses and mass ratios are difficult to estimate for the 
VLM objects in the AstraLux survey. The dominant dif- 
hculty stems from the fact that the ages of the stars are 
only very loosely constrained, which means that the tran¬ 
sition from a star’s brightness to its mass will be accom¬ 
panied by very large uncertainties. Furthermore, the evo¬ 
lutionary and atmospheric mod els (|Hauschildt et all 1 19991 : 
iBaraffe et al.l[200^ : lAllardl[20l4 l that are used to make such 
a transition are themselves uncertain, since they have not 
been calibrated against observations for large sections of the 
parameter space. Both of these issues are particularly criti¬ 
cal for masses that approach the BD range. This is why the 
observed mass-ratio distribution should not be too much re¬ 
lied on. Regardless, as already stated, the observed features 
of interest for our purposes are present in the data at a sta¬ 
tistically significant level. 

In the DPS model, mass ratios for BD-like objects prefer 
values close to unity, as observed for BDs in the Galactic 
held. The DPS model facilitates this through biased-pairing 
fSection 13.21) . Fig. |4] demonstrates how combining this with 
star-like M-dwarfs leads to a distribution that increases with 
increasing q and hattens beyond q « 0.5, for binaries in the 
semi-major axis range log]^g(a/AC/) = —0.4 to 1.95. 


M-dwarfs as tracers of star and BD formation 7 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

primary mass m, / 


□ 

BD desert?- 


X AstraLux late-M (SC sample) 
AstraLux early-M (CS sample) 
* MinMS sample 
- □ HST/NICMO^S 


Figure 5. Primary vs. secondary component masses in the 
AstraLux jjanson et ^ I2OI2I . [20141) and MinMS M-dwarf data 
dWard-Duong et al.l 120151). as well as the HST/NICMOS VLM 
star data (squares. IPieterich et all I2OI2I I. Color-coding as in 
Figs. m and |3] The DPS model prediction is depicted as grey 
dots. The density of dots is a measure of the expected numbers. 
VLM M-dwarfs and BD binaries have a tendency toward unit 
mass-ratio. The upper solid triangle contains binaries from the 
star-like branch, and the lower dashed triangle contains binaries 
from the BD-like branch of the IMF. Where the triangles overlap, 
binaries from both branches can be found. The horizontal dotted 
line is the hydrogen burning mass-limit which roughly marks the 
companion detection limit for the early M-dwarf surveys. 


4.3 A BD desert in the data? 


If BDs form in the same way as stars the natural expec¬ 
tation is that BDs pair in the same way with stars as 
stars do among each other. And if the IMF is continu¬ 
ous over the hydrogen-burning mass li mit BDs should be 
the m ost abundant companions to stars llKrouna fc Bouvierl 
I 2 OO 3 II . However, studies of VLM objects have demon¬ 
strated that this is not the c ase, i.e. a BD desert is ap¬ 
parent (i McCarthy et al.l l2003l : iGrether fc LineweaverllioOfil : 
Dietench et_al.l l2012fl . The SPH plus Nbody approach of 


Li et al.l ~ 2015l l to study the BD-like formation mode, i.e. 


disc fragmentati on, is in agreement with the observed BD 
desert. Although IGrether fc Lineweaverj (I2OO6I ') suggest that 
a BD desert is a natural con sequence of a u n iversa l com¬ 
panion mass function (GMF). IPieterich et al.l (120121 ') show, 
using their HST/NICMOS sample, that VLM stars have 
a tendency toward unit mass ratio, i.e. BD companions are 
rare if not absent. This suggests that deviations from a uni¬ 
versal CMF likely exist, at least for the BD regime (see their 
§6.3 and §6.4 for a discussion). 

Being based on random pairing for star-like objects and 
biased-pairing for BD-likes, the DPS model suggests a simi¬ 
lar BD desert to be present. Fig.l^shows that the DPS model 
prediction is overall consistent with the AstraLux and 
MinMS data. The more the primary component mass ap¬ 
proaches the hydrogen-burning limit from above, the closer 
the m ass-ratio gets to unity . This is consistent with the find¬ 
ing of lPieterich et aP (I2OI2I L As the primary mass increases. 
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Figure 6. DPS model bina ry fractions and th ose observed for 
wide binaries in the study of iDhital et alj ||2010| . estimated from 
their fig. 16). The semi-major axis range covers logjQ(a/A17) = 3 
to 5. Beyond spectral type M5 the binary fraction is essentially 
zero, both in the study and in the DPS model. Observed binary 
fractions are lower limits. 

the DPS model suggests that systems with low mass ratios 
should exist, but they are not seen in the late AstraLux ob¬ 
servations. This might reflect the difficulty to detect sys¬ 
tems where the brightness of the primary outshines a po¬ 
tentially present lower mass star or BD. Thus, the DPS 
model predicts that as instrument sensitivity increases, the 
observed low mass-ratio systems should become apparent for 
late M-dwarfs, increasing the observed binary fraction and 
bringing it closer to that predicted by the DPS model. Such 
systems have been fo und in the HST/NICMOS sample 
llDieterich et aLll2012l . squares in Fig. O- 

Whether a similar BD desert exists for this late M-dwarf 
sample is thus difficult to assess with the AstraLux data. If 
a desert were to exist, it would contradict the suggestion that 
the semi-major axis distribution constructed from the very 
same data (Figs. [1] and |3} is a tracer of a common, canonical 
formation mode for all stars and BDs. Detailed observational 
investigations of the BD desert for VLM stars in conjunction 
with their separation distributions will thus provide impor¬ 
tant insights into the issue of the (non-)continuity of star 
formation. 


4.4 Complementary M-dwarf studies 

The SLoWPoKES survey (iDhital et al.l[201^ investigates 
binarity of low-mass, wide common-proper motion binaries 
in a catalog from the Sloan Digital Sky Survey (SDSS). In 
Fig.[^we compare their observationally deduced wide binary 
fractions, a = 10^ — 10® AU, with the DPS model. A similar 
declining trend towards later spectral types is seen for both 
the DPS model and the observations. In the DPS model 
context this is due to the easier break-up of binaries with 
lower mass primaries in their birth clusters since their bind¬ 
ing energy decreases and the uneven mass-ranges covered 
in e ach spectral subclass (shri nking towards later M-dwarfs, 
e.g. iBaraffe fc Chabrieilll99^ ). The DPS model appears to 
overestimate the published binary fraction fo r the earliest 
M-dw arfs by about an order of magnitude. iDhital et al.l 
ll2010l) state that their wide binary fraction is likely a lower 
limit because observational biases and incompleteness play 
a significant role. The excess of model binaries might be re¬ 
lated to this. If not, then this might indicate a short-coming 


of the DPS model. We note, however, that for M3 and later 
spectral types, the DPS model and observed fraction agree 
better. 

At th e small sem i -majo r axis end the spectroscopic 
study by IClark et al.l (l2012l) finds a binary fraction of 
2.9!tQ 8 per cent for cool M-dwarfs binaries from the SDSS 
having log]^g(a/AC/) < —0.4 assuming a uniform prior semi¬ 
major axis dis tribution (unli ke the rising semi-major axis 
distribution of iKrou 3 Il995l) . This value is in agreement 
with the DPS model fracti on of 2.4 pe r cent in the range 
— 1 < log^Q(a/A?7) < —0.4. [Oark et akl ll2012h additionally 
show that the spectroscopic binary fraction is a continuously 
increasing function of primary mass, from BDs to massive 
0-type stars, which is an additional test of the DPS model. 
Since the present work is on M-dwarfs, we investigate this 
dependence in a forthcoming contribution. 

5 DISCUSSION 

5.1 A case for separate populations 

The sole fact that the most straight-forward strategy to com¬ 
bine star-like and BD-like objects is able to reproduce the 
observations is remarkable. If one had thought of the impli¬ 
cations of a separate BD-like population on the late M-dwarf 
separation distribution before the AstraLux data became 
available, the DPS model described in Section [3] using the 
parameters in Tab. [T] would have been devised in exactly 
the same way. No attempt at trying to merely reproduce the 
observations has been made here. Having said this, the DPS 
model thus post-hoc predicts both the AstraLux late M- 
dwarf binary fraction and semi-major axis distribution, the 
possibly different distribution amplitudes between the early 
and late M-dwarf data, and the MinMS data for early M- 
dwarfs. 

The match probability of « 30% is already good. It 
might increase by relaxing the assumption that the whole 
BD-like population, which includes late M-dwarfs from the 
BD-like branch, follows the BD separation distribution in 
the Galactic field. A proper initial distribution for BD-like 
objects in embedded star clusters which has evolved along¬ 
side the star-like objects through dynamical processing is 
not yet available. But since VLM M-dwarf binaries from the 
BD-like branch are on average more strongly bound than bi¬ 
naries with a BD primary, it is to be expected that in reality 
M-dwarf binaries from the BD branch with a larger separa¬ 
tion can better survive dynamical processing in embedded 
clusters. Their separation distribution in the Galactic field 
will thus likely extend to somewhat larger semi-major axes 
than assumed here, and the amplitude of the peak will in¬ 
crease as well. This would, in turn, improve the consistency 
between the DPS model and the observations. 

The DPS model requires an overlap of the star-like and 
BD-like branches of the IMF (Fig. [2}. For the present mod¬ 
elling the BD-like branch is assumed to extend to Q.2Mq, 
as empirically determined for IC 348 (Sec. 13.211 . which is 
consistent with the mass spectrum arising for objects form- 
ing through disc fragmentation i n the SPH computations of 
Stamatellos fc WhitworthI (l2009h . If the cut-off mass for the 
BD-like branch is lowered to, say, « 0.12Mq (the best-fit 
value for Tau) the amplitude of the peak in Fig.|3]is lowered 


DPS model i t 
Dhital et al. (2010, est.) Q 
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as well. This is due to the smaller overlap of the BD-like 
branch with the observed range of primary masses in the 
late AstraLux data which extends down to « O.IMq only. 
The match probability then changes to 11.6 per cent, thus 
discarding the DPS model is possible at less than 2a con- 
hdence. Staying within the DPS model context, the lower 
match probability using a lower cut-off mass as determined 
for Tau might simply imply that Tau-like binary formation 
is not dominant in contributing to the held, as is addition¬ 
ally evidenced by its super-held binary fraction for long- 
period binaries llKohler fc Leinertlll99a : iKraus et al.ll201h . 
see Sec. E2I. Though this cut-off mass is probably not sig- 
nihcantly larger, we note that the uncertainty associated 
with this empirically-determined parameter implies that a 
somewhat larger value is also consistent with the data (E21), 
and this could improve the agreement with the DPS model. 

Uncertainties in the binary fraction of Galactic held 
BDs and in the peak and width of the BD separation distri¬ 
bution are not considered here. In principle, with the avail¬ 
able observations we have the freedom to vary these param¬ 
eters within the uncertainties to further improve agreement 
with the data. However, we refrain here from arbitrarily 
varying these parameters to obtain a better match since the 
reasonable agreement with this simplest DPS model speaks 
for itself. Instead the issue will be revisited once constraints 
for an initial BD-like population become available. 

5.2 On the meaning of “Universality” 

The notion of universal binary formation is that star for¬ 
mation leads to invariant formal distribution functions due 
to physical processes like energy and angular momentum 
conservation and the chemistry of molecular clouds, all of 
which are the same everywhere, except perhaps in very in¬ 
tense star bursts. This implies an environment-independence 
of binary formation (Sec. 15.2.HI , potentially similar to the 
idea that universal star formation leads to an invariant for¬ 
mal distribution of stellar masses, the IMF. These distribu¬ 
tion functions are parent distribution functions, from which 
a particular case is discretised, or rendered. Universal bi¬ 
nary formation does not preclude, however, a dependence 
on the primary component mass, e.g. different birth bi¬ 
nary distributions for BDs and stars and in different mass 
ranges (0.06 — I.SMq, 1.5 — IOMq, > lOM©) are possi¬ 
ble (Sec. 15. 2. 211 . The universal distribution functions can 
be primary mass-dependent by, e.g., adding further forma¬ 
tion channels, like circumstellar disc fragmentation for BDs, 
which do not change the underlying universal physical pro¬ 
cesses. We need such functions to initialise, for example, 
Nbody models in order to study how young and old clus¬ 
ters evolve into the held and associations. 


5.2.1 Environment (in-)dependence? 


An environment-independence of binary formation has often 
been suggeste d but the observational data don’t appear to 


be conclusive ([KrouDa[[201l[: Krouna & Petr-Gotzens[ 

2011 

Marks & KrouDa[[2012l: [King et al.[[2012l: [Marks et al.l 

2014 

Parker 2014: Leigh et al.[[2014f). 


[King et al.[ (2012f) found the orbital-separation distri- 


buttons in seven young star formation regions to be statis¬ 
tically indistinguishable, consisting mostly of hard binaries 


given the regions’ presently observed conditions. By implic¬ 
itly assuming that the regions were not denser in the past 
they concluded that the observed distributions resemble the 
ones at birth. On the other hand, iMarks fc Kroup^ (l2012ll 
demonstrated that these same regions are consistent with 
dynamicall y processed u niversal birth distributions for late- 
type star s jKroupalll995ll if t hey were signihcantly denser in 
the past. [Marks et ah 1 2014ll later re-visited this issue, and 
argued that the two competing scenarios are statistically in¬ 
distinguishable, given the low number of observed binaries 
in these regions. It follows that solutions to the observations 
are degenerate. In order to settle the issue of whether or 
not star formation in these regions was the same with large 
initial binary fractions, we need to know if these regions 
were significantly denser in the past. With this in mind, 
[Marks et al.l (l2014h offered several arguments to help con¬ 
strain the initial cluster densities, and concluded that such 
a scenario is indeed plausible. 

[Parked (l2014l f offered an intriguing suggestion to break 
this density degeneracy. It involves measuring a star form¬ 
ing region’s degree of substructure, and asking if it can 
res ult from dense initial conditions. The author assumes 
the[K rou ^ (Il995h primary-mass independent initial separa¬ 
tion distribution for late-type stars that is evolved in dense, 
initially substructured and subvirial/collapsing clusters us¬ 
ing Nbody models. The simulations produce clusters that 
are too centrally concentrated, compared to actual observa¬ 
tions of star formin g region s (exce pt maybe the Orion Neb¬ 
ula Cluster, ONCf . [Parked ll2014l '> Hnds initial density esti¬ 
mates from a comparison to the region’s observed structure 
that a re smaller than those constrained bv iMarks fc Kroupal 
1 I 2 OI 2 II from a comparison to the region’s observed binary 
separation distributions. Unfortunately, the author does not 
compare the resulting separation distribution in his interme¬ 
diate density computations - which reasonably match the 
present-day stru cture - to the observed distributions. Inter¬ 
estingly, fig. 1 in [Parked l[2014h suggests that these computa¬ 
tions produce separation distributions that lie between those 
observed in Tau or CrA (i.e. an exce ss of long perio d bina¬ 
ries compared to the field, as in the [Krou 3 ([l995fl distri¬ 
bution) and the ONC (i.e. depleted in long-period binaries). 
This is where the observed separation di stributions in the 
rema i ning investigated r egions lie as well ([Marks fc KrouT)a[ 
[ 2 OI 2 I : [Marks et al.[[20i3l . If the bin ary pop ulations in the 
intermediate density computations of [Parken ([2014['> were in¬ 
deed to resemble the observed ones, as we here suggest, they 
seem to support the hypothesis of environment-independent 
birth distributions for late-type stars. This is because the 
computations then simultaneously reproduce the observed 
present-day structure and binary populations . Thi s raises 
the qu estion: Why do [Marks fc Kroupa[ ([2012[ ~) and [Parker! 
([2014) find different initial cluster densities? One possi¬ 
ble contributing factor is the different initial cluster setups 
(spherical and virialize d vs. substructured and collapsing, 
respectively). However, [Parker et al.l l[201lll conclude that 
this should not make any difference as far as the dynamical 
processing of the binary population is concerned. 

[Leigh et al.[ ([2014[ f confirm a similar density degener¬ 
acy of the primordial binary population for globular c l usters 
(GCs). Using MOGCA computations l[Giersz et al.[ [2013fl 
over a Hubble time they show that solutions to reproducing 
simultaneously the rather low binary fractions and an anti- 
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correlation between the binary fractions and masses of GCs 
inside the half-mass radius are degenerate in terms of the ini¬ 
tial GC densities and initial binary fractions, quite similar to 
the density degeneracy for young star forming regions. How¬ 
ever, they break this degeneracy by demonstrating that only 
dense initial configurations, which match densities observed 
for potential GG progenitors, in combination with large ini¬ 
tial binary fractions account for a similar anti-correlation 
seen outside the half-mass radius in Galactic GGs. Thus, 
the observations for GGs are consistent with t he universal¬ 
ity hypothesis (n ot precluding different origins, iLeigh et al.l 
|2012|. 12011 l2014l'l . 


5.2.2 Mass-independent or field-like binary formation? 


The present work suggests that primary-mass independent 
birth distribution functions for late-type binaries plus dy¬ 
namical processing continues to be a valid working hypoth¬ 
esis for M-dwarf binaries. Whether the birth distribution 
functions can be primary-mass dependent for late-type bi¬ 
naries or not is being debated. _ _ 

In a recent contribution, IParker fc MeveJ (l2014l ~) find 
that their Nbody computation with fractal initial conditions 
and peak densities around IOOOMq pc“® do not reproduce 
observat i ons o f binaries in the Galactic field, when using the 
iKroupal lll995ll distribution as input and assuming BDs as 
a continuation of stars with the same universal init i al con - 
ditions. This has been done before in iKrouna et al.l (l200l) . 
who show that this leads to distributions inconsistent with 
observations, and in particular too many star-BD binaries. 
iKroupa et al.l ll200ll concluded that BDs need to be treated 

with their own pairing rules. _ 

_ More precisely, in the computations of[Parter|&^^e^3 

(l2014li . the orbital separation distributions after 10 Myr of 
dynamical processing do not match the Galactic field sepa¬ 
ration distributions for 


(i) brown dwarfs llThies fc Kroup^l2007l '). 

(ii) M-dwarf b i naries in the AstraLux survey 

djanson et al. I l20lll2014lb __ 

(iii) solar-type binaries dRaghavan et al.ll2010ll and _ 

(iv ) A-type binaries in the VAST survey I De Rosa et al.l 
l2014l) . 


Based on this, IParker fc Meveil d2014l l suggest that the ini¬ 
tial distributions are primary-mass dependent. Thus, com¬ 
putations with the observed distributions as input to the 
otherwise same cluster provided a better match to the obser¬ 
vations, although not perfectly due to the dynamical break¬ 
up of wide binaries among the G- and A-type population. 
They conclude that the field binary populations are indica¬ 
tive of the star formation process in clusters. 

The results presented in this paper are consistent with a 
different interpretation, however, which have the advantage 
of being predictive. In the following we address the apparent 
disagreement with the four populations above: 


(i) a meaningful comparison of the dynamically processed 
universal semi-major axis distribution for late-type stars 
and the BD binary population in the Galactic field re¬ 
lies on the a ssumption that star formation is continuous 
and that the l& oupal (Il995ll distribution is valid for BDs 
dParker fc Mever 20141 '). If this were not true, the compari- 



Figure 7. Comparison of orbital period distributions for solar - 
type stars (histogram with Poisson errors. IPaghavan et ahlboKt) 
and G-dwarfs in the DPS model (thick crosses). The same param¬ 
eters (Tab. [TJ used to construct the M-dwarf DPS model (Fig. (Sjl 
are employed. 


son would not be meaningful. As d emonstrated here, re quir¬ 
ing a separate BD population as in iKroupa et al.l ll2003l) and 
iThies fc KrouPc] (l20()7l . l2008l ) yields agreement with the ob¬ 
served BD binary semi -major axis dis tribution by construc¬ 
tion. Importantly, the iKroupal lll995h distribution was not 
developed to match the observed properties of binary BDs, 
but was instead developed for (originally only) solar-type 
stars after dynamical processing. 

(ii) as shown in this contribution, t he unive r sality hy¬ 
pothesis for late-type stars following the iKrou ^ (Il99^ or¬ 
bital parameter distributions, combined with dynamical pro¬ 
cessing and a separate BD-like population, can successfully 
reproduce the late AstraLux and early MinMS M-dwarf 
data. 

(iii) extracting the orbital period distribution for G- 
dwarfs in the Galactic field from the DPS model adopting 
the same parameters used to extract the M-dwarfs (Tab.[T} 
shows agreement between the observed pe riod distribution 
of so lar-type binaries in the Galactic field (|Raghavan et al.l 
l2010l. including orbits in hierarchical multiples) and the 
DP S model (Fig. [3). This is not surprising, though, since 
the iKrouoal | l99^) model was designed to match the solar- 
type data of iDuauennov fc Mavoil (ll99ll L and their sepa¬ 
ration distributi on of orbits is indistin guishable from the 
one obtained by iR.aghavan et al.l (l201f)l). On the or i gin o f 
the different conclusions reached bv IParker fc Me^ J2014h 
we can only speculate. We note that Kroupa ( 19951 ) ex- 
plicitly demonstrated that only certain combinations of ini¬ 
tial cluster masses and radii reproduce the observed pre- 
mainsequence and field binary populations simultaneously. 
This yields constraints on the physical conditions in corre¬ 
lated star-forming events (i.e. embedded clusters). 

(iv) while A-type binaries in the VAST survey show a 
narro w distribution peaked around 370 au dPe Rosa et al.l 
l2ni4h . spectroscopic A-type binaries a r e also known to ex- 
ist at shorter separations dAbtl Il965l : iGarrier et al.l l2002l : 
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ICarauillat fc Prieuill200^ . Even a doubl e-peaked distribu¬ 
tion for A-type binaries is suggested dPuchene fc Kraua 
l201.'j ) . However, spe c trosco pic data is needed to comple¬ 
ment IPe Rosa et alJ (l2014l i’s data for a complete picture. 
We will address the VAST data in future work. However, 
the (initial) distributions for O a nd B-type bina ries appear 
to be different from that given in iKrou 3 ll 19951 1. This may 
be indicative of a different formation channel for O- and 
B-type stars, e.g. through the star-fo rmation process at th e 
density centres of proto-clusters fe.g. lBonnell &: Batell2005ll . 
analogous in some ways to a separate BP formation chan¬ 
nel. In this case, the PPS model should be extended to ac¬ 
count fo£_observations of binaries with hi gh-mass p ri marie s 
(e.g. ISana et ^ l2012fl . as developed by loh et al.l ll2015h . 
However, one would naively expect any change in the star 
formation process to be mapped onto the stellar IMF (cf. 
iBaneriee fc Kroupall2012ll . 


The take-away message is that both a field-like for¬ 
mation scenario (“what-you-see-is-what-you-get”) and the 
PPS model with a separate BP-like population explain the 
Galactic field observations. The latter formulation, assuming 
one birth binary population distribution function for low- 
mass stars (< I. 5 M 0 ) plus dynamical processing to be the 
origin of all late-type stellar binaries and adding a sepa¬ 
rate BP-like population (and perhaps a separate OB pop¬ 
ulation) explains the data. This scenario is not only con¬ 
sistent with, but actually explains the large range of esti¬ 
mated global__biMn;_^;actions, from < 10 per cent in oj Cen- 


tauri telson et ah 1995h to sub-field fractions in other GCs 
llMilone ef^ah 2012fl to > 90 per cent in young s tar forming 


regions I Kohler fc Leiner3 Il998l : TPuchenel Il999l ~l. Observa¬ 


tions of separation distributi ons in the latter objects might 


iMarks & KrouDall2012l:lMarks et al. 

2 OI 4 I'). In addition, very 

young regions like Tau and CrA ( 

Kohler & Leinerd 119981: 

iKohler et alJl2008l:lKraus et al. 

l201lh as well as oroto-stellar 

binaries (Connellev et al.l 2006 

) exhibit a binary excess at 


long periods that do not support a field-like formation sce¬ 
nario. These populations all have a significant soft compo¬ 
nent, which changes quickly through dynamical processing 
in intermediate to dense environments. Such observations re¬ 
sult naturally from different degrees of dynamical processing 

of the underlying binary populations. _ 

Thus, the dynamical modification of the lKrou ^ (Il995ll 
distribution combined with a separate BP-like population is 
not only consistent with the data, if the initial cluster density 
is high, but it also allows binary populatio ns to be predicted 
in th e Milky Way and other galaxies llMarks 8z Kroupa 
l201lh . The binary excess in the younge st and sparsest star 
forming regions is accounted for by the iKroupaf (Il995h dis¬ 
tribution used in the PPS model, and its dynamical mod¬ 
ification yields the bi nary populations in star clusters and 
star f orming regions llMarks etliL 201l|_ Marks fc Kroup^ 


l2Olj;jMnjk£ef0^j2Olj: iLeigh et al.! 12014 1 and in galaxies 


( Marks fc Kroup l201ll. and this contribution). 

We caution that taking observed distributions and us¬ 
ing them as input for computations of clusters which re¬ 
semble the presen^ observed state as initial conditions (or 
even a moderateljO dense configuration) will in most cases 


not change the observed binary distributions significantly, 
unless the observed cluster is extremely young (^ IMyr) 
such that the binary population did not have the time yet 
to adapt to its host cluster properly. This is because dy¬ 
namical binary processing occurs rapidly if the population 
had been initially binary dominated and will rather quickly 
find an equilibrium with its host environment (the changes 
happen on a crossing-time scale, iMarks et al.ll201ll ~) before 
the long-term, two-body relaxation driven evolution domi¬ 
nates further processing llLeigh et al.ll201^ . Thus, such an 
equilibrium state is what one will most likely observe. 


5.3 Predictions 


Whether BPs and stars are indeed separate populations 
with their own initial mass and binary distribution func¬ 
tions can be constrained with future surveys of Galactic 
field M-dwarf binaries. The presence of a separate BP-like 
subpopulation unmasks itself through a peak located close 
to the field BP separation distribution which resides on- 
top of a wider distribution, as a result of the superposi¬ 
tion of the star-like and the BP-like M-dwarfs (Fig. [S] left 
panel). The PPS model predicts that M-dwarf binaries in 
the field exist in significant numbers outsi de the sensitivity 
range of the AstraLux s urvey (e.g. as in iFischer fc Marcvl 
ll992l : |Pelfosse et al.ll2004h . i.e. the separation distribution is 
wider than reported by the AstraLux data, bo th for early 
djanson et ahUioi^ and late (Ijanson et al.ll2 0l4l M-dwarfs. 
The recently published MinMS data 1 Ward-Puong et ^ 
I 2 OI 5 II supports such a notion for the early M-dwarfs. If 
true, the late M-dwarf data roughly reproduce the peak 
of the PPS model distribution, whereas the observed early 
M-dwarf data missed it slightly due to the sensitivity con¬ 
straints. 

Here the BP-like branch of the IMF extends up to 
0.2Mq. All M-dwarf binaries with < 0.2 Mq in the « 1 — 
10 AU range are thus apriori BP-like binary candidates. By 
comparing these objects among each other, observers might 
be able to find further signatures for a BP-like formation 
mode, aside from the peak in the separation distribution. 
These signatures might be hidden, e.g., in other orbital pa¬ 
rameters such as their eccentricity. If the true BP-like M- 
dwarf binaries follow the Galactic field BP mass-ratio distri¬ 
bution they preferentially have a mas s-ratio close to unity. 
The analysis of disc fragmentation in iLi et all ll2015ll pro¬ 
vides further clues. 

We expect a peak superposed on a wider distribution 
in M-dwarf binary populations which have primary com¬ 
ponents with masses below 0.2Mq. The peak is not ex- 
pected if such bina ries are absent, as is e.g. the case for 
Ijanson et al. I ll2012li’ s e arly M -dwarfs and the MinMS sam¬ 
ple I Ward-Puong et al.ll2015l i. Their sensitivity limits do not 
include the separation range where the peak is expected and 
prevents us from testing this scenario. Should a peak be seen 
for early M-dwarfs in later surveys as well, the maximum 
mass of BP-like objects needs to be shifted to larger masses. 
This provides thus a means to constrain observationally how 


below about the maximum of the input semi-major axis distribu- 
® I.e. a configuration which does not place the hard-soft boundary tion. 
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deep the separate BD-like population, if it exists, penetrates 
into the star-like regime. 

Our random pairing scenario of binary components 
for star-like M-dwarfs predicts that low mass-ratio late 
M-dwarfs will be detected at larger primary masses and 
that these will be lacking BD companions (Fig. [Sjl, i.e. ex¬ 
tend the BD desert. These are currently invisible due to 
AstraLux instrumental sensitivity constraints. 


6 SUMMARY 

It has been demonstrated how combining a dynamically pro¬ 
cessed, initially binary-dominated universal star-like popula¬ 
tion for late-type stars (< 1.5M©) with a separate BD-like 
population forming through circumstellar disc fragmenta¬ 
tion, which produces some VLM stars, leads to a prediction 
for the semi-major axis distribution of late M-dwarfs in the 
Galactic field. Upon adding the two formation channels . Dy¬ 
namical Population Synthesis (iMarks Sz KrouI)3l2011^ pro¬ 
duces a narrowly peaked distribution that appears on-top 
of a wider distribution. The simplest DPS model resembles 
the late M-dwarf d ata obtained in the AstraLux survey 
djanson et al.|[20l3 j and the null hypothesis that the DPS 
model provides a parent distribution for the observations 
is confirmed to within Icr certainty. It has been discussed 
that agreement between DPS model and observation could 
improve once an estimate for an initial BD-like binary pop¬ 
ulation becomes available and once it is allowed to partici¬ 
pate in dynamical processing alongside the star-like objects 
in the DPS model. As a byproduct our results demonstrate 
that the hypothesis of primary-mass independent period, or 
semi-major axis distribution functions accounts for all late- 
type populatio ns in the Galactic field. The opposite conclu¬ 
sion reached bv IParker fc Meveil ll2014ll does thus not follow 
unambigiously from the data (see Section 15.2.21 for a discus¬ 
sion) . 

The DPS model predicts that binaries outside the sen¬ 
sitivity range of the AstraLux survey will be detected. A 
peaked semi-major axis distribution superposed on a wider 
distribution as for the late M-dwarfs is not expected for early 
M-dw arfs (as those in ijanson et al.ll2012l : IWard-Duong et al.l 
l2015h . Once surveys for early M-dwarfs become sensitive 
to smaller separations than covered by the AstraLux and 
MinMS data this notion will provide a means to measure 
how deeply the BD-like branch penetrates into the stellar 
regime. The DPS model separation distribution agrees ex¬ 
cellently with the early M-dwarfs in the MinMS sample 
but less well with the early AstraLux data. It is pointed 
out that the ob served separation d istrib utions for early M- 
dwarf s found bv ijanson et al.l (l2012tl and IWard-Duong et al.l 
(l2015h . respectively, are potentially in tension fSec. 14.iT) . 

Given the surveys of nearby stars, a BD desert is ex¬ 
pected to be present in the AstraLux data. No BD compan¬ 
ion is seen in the data above a primary mass of 0.2M© but 
the decreasing detection probability of a low mass-ratio sys¬ 
tem with increasing primary mass does not yet allow for any 
firm conclusions. If a BD desert in the data were confirmed 
it would contradict the continuous star formation scenario 
inferred from the narrow observed M-dwarf separation dis¬ 
tributions. This is because it is unclear why stars and BDs 
should pair differently, avoiding each other while forming in 


the same way (for a discussion on the observational reali ty of 
the BD desert see §6.3 and §6.4 in iDieterich et aLll2012ll . In¬ 
vestigating the BD desert for VLM stars together with their 
separation distribution in a single sample could thus lead 
us closer to answering the question of the (non-)continuity 
of star formation. 

To summarize, the available M-dwarf data do not pro¬ 
vide unambigious evidence for star formation that is con¬ 
tinuous over and above the hydrogen-burning mass limit. 
This is because, as illustrated in this paper, the data are 
also consistent with the here devised DPS model which 
assumes separate formation modes for stars and BDs and 
universally valid birth separation distributions for late-type 
stars. Specifically, the DPS model is predictive and able 
to reproduce both the MinMS and the AstraLux (with 
poorer agreement) data, as well as the data for other 
binary po pulations observed in a diverse range of envi¬ 
ronment s ([Mar ks fc Kroupall201l] . l2012l : iMarks et al.ll2of3 : 
iLeigh et al.ll2014l L 
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